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Abstract 

To obtain the explicit form of evolution operator in the Tavis-Cummings model we 
must calculate the term e -**9(S'+®a+s , _®a t ) explicitly which is very hard. In this paper we 
try to make the quantum matrix A = S+ &> a + S- (g> at diagonal to calculate e~ ltaA and, 
moreover, to know a deep structure of the model. 

For the case of one, two and three atoms we give such a diagonalization which is first 
nontrivial examples as far as we know, and reproduce the calculations of e~ ttgA given in 
quant-ph/0404034, We also give a hint to an application to a noncommutative differential 
geometry. 

However, a quantum diagonalization is not unique and is affected by some ambiguity 

arising from the noncommutativity of operators in quantum physics. 
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Our method may open a new point of view in Mathematical Physics or Quantum 
Physics. 

1 Introduction 

The purpose of this paper is to give a new insight to the Tavis-Cummings model ( lj) and to 
obtain the explicit form of evolution operator by the new method in the case of some atoms. 

This model is a very important one in Quantum Optics and (maybe even) in Mathematical 
Physics, and has been studied widely, see [2] as general textbooks in quantum optics. 

We are studying a quantum computation and therefore want to study the model from this 
point of view, namely the quantum computation based on atoms of laser-cooled and trapped 
linearly in a cavity. We must in this model construct the controlled NOT gate or other controlled 
unitary gates to perform the quantum computation, see as a general introduction to this 
subject. 

For that aim we need the explicit form of evolution operator of the model in the case of (at 
least) one, two and three atoms. As to the model of one atom or two atoms it is more or less 
known (see [4J), while as to the case of three atoms it was given by jH]. However, the method 
is not clear enough in a mathematical sense 1 . 

In this paper we present a quantum diagonalization method which is a quantum version of 
classical diagonalization and obtain the explicit form of evolution operator obtained in j^J. 

However, the quantum diagonalization is not unique and is affected by some ambiguity 
due to the noncommutativity of operators in Quantum Physics. This may be related to the 
so-called operator ordering problem, see for example P on this topics. 

The Quantum Diagonalization Method is completely new and may be applied to a noncom- 
mutative differential geometry (for example, the noncommutative chiral models) because we 
can construct (quantum) unitary matrices explicitly. However, this is beyond our scope of the 
paper. 

1 We used Mathematica in the process of calculation 
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2 Tavis— Cummings Model and Evolution Operator 

We make a review of [H] within our necessity. The Tavis-Cummings model (with n-atoms) 
that we will treat in this paper can be written as follows (we set % — 1 for simplicity). 

A n n 

H = uj\ l a* a + - ® 1 + 9j2{ a i +) ® a + ^ ® flt ) > W 

^ i=i i=i 

where u is the frequency of radiation field, A the energy difference of two level atoms, a and 
a) are annihilation and creation operators of the field, and g a coupling constant, and L = 2 n . 
Here u\ + \ u\ ' and af^ are given as 

o\ s) = h ® ■ • • ® I2 ® o- a ® 1 2 ® • • • ® I2 (* - position) e M(L, C) (2) 
where s is +, — and 3 respectively and 

H;;)' '-=(;;)• Mi -J- i2 M°)' <3) 

Here let us rewrite the hamiltonian ([Tjl. If we set 

n n 1 n 

s+ = £4 +) , s- = Y,°t\ Sz = ~Y,°?\ (4) 

i=l i=l / i=l 

then (JJJ can be written as 

H = ool L <g> a f o + AS3 <g> 1 + g (s+ ® a + S_ ® a 1 ") = # + V, (5) 



which is very clear. We note that {S + , SL, 5*3} satisfy the su(2)-relation 

[5 3 , 5+] = S+, [S 3 , SL] = -S-, [S+, 5_] = 25 3 . (6) 

However, the representation p defined by p(o+) = S + , p(er_) = S_, p(a 3 /2) = S3 is a reducible 
representation of su(2). 

We would like to solve the Schrodinger equation 

ijU = HU= (H + V) U, (7) 



where U is a unitary operator (called the evolution operator). We can solve this equation by 
using the method of constant variation. The result is well-known to be 



under the resonance condition A = u, where we have dropped the constant unitary operator 
for simplicity. Therefore we have only to calculate the term (JSJ) explicitly, which is however a 
very hard task 2 . In the following we set 

A = S + <g> a + S_ <g> a ] (9) 

for simplicity. We can determine e~ ltgA for n = 1 (one atom case), n = 2 (two atoms case) and 
n = 3 (three atoms case) completely. 



One Atom Case In this case ^4 in (jOJ) is written as 

a 
v a) 

By making use of the relation 

' aa) \ I N + 1 
a^a ) \ N 
with the number operator iV we have 

cos^v^VTT) -,^±1) 

(tgVN) 



(10) 



At 



-itgAi 



-la 



:i2) 



VN+l 



COS 



We obtained the explicit form of solution. However, this form is more or less well-known, see 
for example the second book in [2] . 

Two Atoms Case In this case A in 



A, 



V 



is written 


as 


a 


a 





a) 





a 


a) 





a 


a f 


at 






(13) 



/ 



2 The situation is very similar to that of the paper quant-ph/0312060 in 7 



Our method is to reduce the 4 x 4-matrix A 2 in ()13|) to a 3 x 3-matrix B\ in the following 
to make our calculation easier. For that aim we prepare the following matrix 



/ 



10 



75 
75 



1 



then it is easy to see 



T^AoT 



V 



V2a 
v / 2a t 
V2d 








(14) 



where Si = J + ® a + J_ <S> a' and {J + , J_} are just generators of (spin one) irreducible repre- 
sentation of 0. We note that this means a well-known decomposition of spin | ® ■ 



Therefore to calculate e %tgA2 we have only to do e ltgBl . Noting the relation 



/ 



2(2iV + 3) 



\ 



2(2iV + 1) 



B = DBi, 



2(2iV-l) y 



we obtain 



-itgBi 



l + t£±if(AT+l) -i/i(JV+l)a 

1 + 2/(iV) 



2iV+3- 

-m t /i(iV + 1) 



-ih(N)a 



\ Ww+ 3 f(N + 1) -ia^h(N) 1 + ^/(JV - 1) j 



where 



f(N) 



1 + cos (tgJ2(2N + 1 



sin [tg^2{2N + l] 
V2iV+ 1 



Offil. 



(15) 
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Three Atoms Case In this case A in is written as 



a a 
at 



a 



a 
a a 
a a 



at at a 
atOOOOaaO 
at at a 
at a f a 
at at at 



(16) 



We would like to look for the explicit form of solution like f!12D or p5j) . If we set 



i 

-1 
V2 
















1 

-1 

V2 



10 



1 

76 

1 

75 









73 
73 

% 

i 

75 

















1 



1 

V3 



- 73 



n/6 





1 



then it is not difficult to see 

^ a 
at 



a 
at 





V^at 





\^3a 


2at 






2a 


\/3at 










B 



1/2 



B 



1/2 



B 



(17) 



3/2 / 



This means a decomposition of spin = \ ®\ ®\. Therefore we have only to calculate 

e -it g B 3/2 ^ which i s however not easy. The result is 



e -itgB 3/2 

( f 2 (N + 2) -VSiF 1 (N + 2)a 2VSh 1 (N + 2)a 2 -QiH (N + 2)a 3 \ 

-V3iF 1 (N + l) a t fi(N + l) -2ii?i(JV + l)a 2V3h(N + l)a 2 

2 v / 3/ii(iV)(a t ) 2 -2ii/ 1 (Ar)a t / (iV) -V5iF (N)a 

\ -QiH (N- l)(at) 3 2 v / 3/ii(A^- l)(at)2 -^F (iV - l) a t f-i(N — 1) ^ 



where 



/ 2 (iV) = 
fi(N) = 
fo(N) = 
f-i(N) = 
hi(N) = 

F 1 (N) = 

F (N) = 
Hi(N) = 
H (N) = 



{v + (N)cos(tg^\ + (N)) - v4N)coB(tgy/\-(N))} /(2y/d(N)), 
{w + (N)cos(tgyJ\ + (N)) - w_(N)cos(tgy/\_(N))} /(2yJd(N)), 
{v + (N)caa(tgy/\-(N)) - v-(N)caa(tgy/\ + (N))} /(2y/d(N)), 
{w + (N)cos(tgy/\-(N)) - w-(N)cos(tgy/X + (N))^ /(2y/d(N)), 
|cos(^ v /A + (iV)) - cos(^^/A_(7V))} /{2^d{N)) ) 

^Lsin(^ V / A^iV : )) - -^^ S in(tg^{N))) /(2^d(N)), 



v + (N) 



Mtgy/\-(N)) 



:sm(tgy/\ + (N))\/(2y/d(N)), 



y\ + (N)sm(tg^X + (N)) - ^A_(JV)sin(t^A_(iV))} /(2y/d(N)), 

1 . , r, TTTTx 1 



^/a^aO 



sin(^ V / A + (iV)) - ^_=sin(t^A_(iV)) [ /^y^) 



and 



(18) 



A±(iV) = 5iV ± }Jd(N), v±(N) = -2N - 3 ± }Jd(N), w±(N) = 2N - 3 ± ^Jd(N), 
d(N) = 16N 2 + 9. 

This form is very complicated. We note that to calculate e~ lt9B:i /' 2 we used Mathematica to 
the fullest. 



3 Quantum Diagonalization Method 



First of all we explain the method which we call a Quantum Diagonalization Method (QDM). 
To calculate e~ it9A for A in © 

A = S + <g> a + S_ <g> a f 

we would like to diagonalize it like A = UDaU^ {Da is a diagonal matrix) if possible. This 
is a well-known classsical procedure. However, in our case it is impossible because we cannot 
determine the eigenvalues by making use of its characteristic equation /(A) = det(Al — A). In 
the quantum case there is no meaning on determinant function. For example, which is correct 

/(A) 

Therefore we have no general method to make A diagonal. However, we have a very skillful 
method for A whose procedure goes like 

Classicalization — > Quantization — > Classicalization. 

The (quantum) matrix A above can be decomposed as 

T* AT = ^2^Bj (a direct sum of quantum matrices of spin j) 
by an orthogonal matrix T 3 like in the preceeding section, where Bj is given by 



A —a 


H 


A 2 - aat 


-at A 




A 2 — ata 



( 



yJ(J-l)la 
(J - l)lat J (J - 2)2a 



(J - 2)2at 







(J (J - k)kaS k ^i + yJ(J - i)ia ] 5 k _ hi 
3 To find T in the general case is not easy 



(J-3)3a 



v /2(J-2)a t 




yJl(J-l)a 
1(J- l)at 



/ 



(19) 
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with J — 2j + 1. In the following we set B = B~ for simplicity. 



(i) Classicalization We replace a — > z and a" 1 " — > 2; in S and set 



C = 
/ 



yl(J-\)\z 
y/(J-l)lz ^(J -2)22 

(J -2)2* J(J-3)3z 



2(J-2)z 
^l(J-l); 



1(J- 1)* 




\J(J — k)kz5 k ,i-i + \]{J - i)iz8 k ^}j . 



We must diagonalize C. The eigenvalues are 



{(J-l)|z|,(J-3)H,...,(J-2i + l)H,...,-(J-3)H,-(J-l)H} 



(20) 



and corresponding orthonormal eigenvectors are 

|(J -2t + 1)\Z\) = (xki-TTj^ 



for % — 1 ~ J 



where x^i are defined as 



2/w 



with y ki (k — 1 ~ J) defined by the recursion relation 



2/1", 



= 1, \j{J -k + l)(k - l)y k -i,i + \/{J ~ k)ky k+hi — (J — 2i + l)y ki . 



For example, 



J-2i + l (J-2i + l) 2 - (J- 1)1 
= 1, 2/2i = / 2/3i = t=_ etc. 



(J - 1)(J -2)1-2 
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We note that the matrix X = (x^) is an (real) orthonormal one, namely X T X = XX T = lj. 
For example, when j = 1 ( J = 3) it is easy to show 



X 



( i i i \ 

J_ n l 

. 1 1_ 1 , 

\ 2 V2 2 / 



If we set 



= 



lfc-1 



(21) 



then W is a unitary matrix (W^W = WW^ = lj) and C is diagonalized by W as 

C = WD C W ] (22) 

where D c is a diagonal matrix consisting of the eigenvalues {(J — 2i + l)\z\ \ % = 1 ~ J}. 
We note that the unitary matrix W is not defined at z = 0, see (J5TJ). 

(ii) Quantization Next we consider a quantization of : namely we want to find a (quan- 
tum) unitary matrix U\ arising from W above. After some trial and errors we set 



then it is not difficult to check 



yjN{N - 1) • • • (N - k + 2) 



(23) 



on the representation space H © Tii © • • • © 'Hj-x where TCk = Vectc{|&), \k + 1), • • • } is the 
subspace of the Fock space H = Ho generated by {a, a*, iV}. We note that U\ is not defined on 
the whole space 4 . 

A comment is in order. Noting 

( a t)V = jy(iV - 1) • • • (N - I + 1), a l (a<) 1 = (N + l)(N + 1 - 1) • • ■ (N + 1), 
aN = (N + l)a, a ] N = (N - l)a f 



1 /7i is a partial isometry on 7i © 7i < 



) Tt in the mathematical terminology 
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for / > 1 we have 



(at)'a' 



(gt) 



1 



N(N -1)---(N -l + l) 
1 



tv 



(gt) 



= a 



(iv + z)(iv + z-i)---(iv + i; 



a'(at)' 



and 



(at) '^y a ' ^ (at) V + 1) . . ■ (at + i - D(n + o "' = 1 on H " 

For example, when j = 1 ( J = 3) we have 

\ 



/ 

2 



1 



1 fl t 

i W 



i 

2 

W 



V 2^(^-1) 



— . 1 (at) 2 . 1 (gt) 2 



(24) 



(25) 



(iii) Classicalization Here we consider a diagonalization of B by U\ above. Some calculation 
leads 

U\BU X =R= (r kt ) (26) 



where 



ik- 



r ki = J2y( J + i)(yi-i,kyu + yi,kVi-i,iW N + 1 - 1 = n, 

1=2 

We see that the matrix R is hermitian and its entries consist of some functions of the number 
operator N, so R is a kind of classical matrix. Therefore we can make R diagonal 5 like 



R = U 2 D R UI 



(27) 



For example, when j = 1 ( J = 3) we have 



R = 



( 



V 



V2 



v / Af+2-y / ]V+T 





v / jV+2-y / ]V+T 
V2 



V2 



3 To obtain Ui explicitly is not easy or almost impossible 



11 



and 



D 



( v /2(2Af+3)+y / ]V+2+v / ]V+T y/2{2N+Z)-y/N+2-y/N+l \ 



2y / 2(27V+3) 



v / jV+2-y / ]V+T 
^^2(2^+3) 



2y/2(2N+3) 



2(2A + 3) 



v / 2 v /2 ( 27V + 3 ) 

v / ]y^+v / An : T 

^2(2^+3) 

^^2(2^+3) 
\ 



V 



2(2A + 3) y 



As a result we finally obtain 



2y / 2(2iV+3) 
V2^2(2N+3) 



v /2(2Af+3)+v / 7V+2+ v ^V+T 
2^/2(2^+3) 



(28) 



This is just the diagonal form of B that we are looking for. We note that all entries ofU = U1II2 
consist of N and a) (not contain a) . From this we have 



e -itgB = UQ-itgDjj]^ 



(29) 



This is a kind of "normal ordered" diagonal expression of the evolution operator. See 
Appendix for another diagonal expression. In the following we give an explicit expression in 
the case of one, two and three atoms. 

Let us list the results. 
One Atom Case For A\ in (fTUj) we have 



A x = UDU^ 



(30) 



where 



U 



/N a - 



D 



I 



V 



-y/WTT J 



Then it is easy to see 



(31) 



e -UgAi = Ue -itgDjj] = the right hand gide of ^ 



f32) 
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Two Atoms Case For B 1 in (JT3j) we have 



where 



U 



^2(2^+3) 



B x = UDU ] 



_y^7V+2 = 

v /2(2AT+3) 





^2(2^+3) 



\ 



,-h2 



V v / iv r V 2 ( 2JV - 1 ) 



fat) 2 



/Ny/2(2N-1) 



and 



Then it is not difficult to see 



^2(2iV + 3) 



ViV-TA/2(2^Ty V ; / 



-,^2(2iV + 3) y 



e -%Bi = Ue -it9Dtf = right hand gide of ^ 



Three Atoms Case For _B 3 / 2 in (fTTj) we have 



where 



U 



B 3/2 = UDU ] 



( \ 

«11 «12 «13 «14 



"21 M22 «23 «24 
«31 M 32 M33 M34 

y m 4 i m 42 M43 M44 y 



where 



u n = j{V2(l-/3 7 ) + V6(/3 + 7 )} 



i{V6(l-/3 7 )-V2(/3 + 7 )} 



M12 



«13 = «12, Ml4 = -Mil, 



^/(l + /?2)(l +7 2)- 
1 



3 The calculation in this case is interesting and hard, so the full details will be given in [5] 
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U21 



\ \^La\l + Pi) + ^La\p - 7) ) 1 



M22 = l(^ Q t(i + / g 7 )-^ Q t( /3 _ 7 )\ 1 

4\VN HfJ v^T 'J ^(1 + ^(1 +7 2)' 

«23 = -U22, U24 = «21, 

7 ( ; = ( a f ) 2 (l-/37) - ; = ( a f ) 2 (/3 + 7) / ■ 
4 I ^ - 1) JN(N - 1) J v /(l + /5 2 )(l+7 2 ) 



U31 



«32 = 7 



at) 2 (l - 7 ) + . vu (at) 2 (/j + 7 ) 



4 \ ^7V(7V - 1) ' 0V(7V - 1) ' J ^/(i + ^)(i+ 7 2)' 

M33 = M32, m 34 = -• u 3 i, 

«4i = )[ -, = (at) 3 (l + /j 7 )- ^ (a t ) 3 (/3- 7 ) 

4 I yfN(N-l)(N-2) ^N{N -1)(N -2) 



x 



M42 



/(l +)9 2)(l + 7 2) 



, , (at) 3 (l + /3 7 ) + v =(a t ) 3 (/3 - 7) 

4 \ ^iV(iV-l)(iV-2) yjN(N-l)(N-2) 

1 

M43 = -M42, = «4i (39) 

and 

H-u-(x-y) fM + u-(x + y) 

f 2b ' 7 2c 

x = 3 ViV + 1 + 6ViV + 2 + 3ViV + 3, y = 3VN + 1 - 2VN + 2 + 3VN + 3, 

b = 2V3 (VnTT - VWT3) , c = Vz (y/WTT - 2VWT2 + VWTs) , (40) 

and 



// = 4]J5(N + 2) + y/16(N + 2) 2 + 9 = 4y/\+(N + 2), 

v = 4^5(iV + 2) - y/l6(N + 2) 2 + 9 = 4y/\-(N + 2) (41) 
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and 



v 



-fi J 



^A + (iV + 2) 



Then we obtain 



e -itgB 3/2 = Ue-itgDtf = the right hand gide of (JJHJ). 



(42) 



(43) 



However, the proof is not easy, see Appendix. 

Last we make one comment. We would like to perform a diagonalization for the case of more 
than three atoms, however it is not easy at the present. One of main difficulties is the step (iii). 
That is, to perform a diagonalization to the (classical hermite) matrix we must determine its 
eigenvalues by solving the characteristic equation and give orthonormal eigenvectors explicitly 
(not abstractly). The characteristic equation is in general algebraic one of degrees more than 
four, which is impossible to solve in an algebraic manner by the famous Galois theory. Even for 
algebraic equations of degrees three and four we must use the Cardano and Ferrarri formulas 
(see for example [9 ) which make hard to determine all orthonormal eigenvectors explicitly. 



4 U(l) Ambiguity 

In this section we discuss a problem of £7(1) ambiguity in the quantum diagonalization method. 

The classical diagonalization C = WDcW^ in ()22jl has the following U(l) "invariance" , 
namely 

C = WD C W ] = WU (U j D c U )U j W^ = (WU )D c (WU o y 



15 



where U is a diagonal matrix defined by 



Un 



G U(J) 



and fj(z,z) is an element in U(l), \fj(z,z)\ 2 = 1. For example, fj(z,z) 
diagonal matrix Dc is invariant under the change of unitary matrix W 
However, this is not kept in the process of quantization. From (|28|) 

B = UDU ] = (UU )(UlDUo)(UU o y = UDU ] 

where Un is a quantum "diagonal" matrix defined by 



ojn. Namely, the 

. wu . 



t 



Un 



/i(a,a f 



f 2 {a,a j < 



V 



and fj(a,a)) is an element satisfying fj(a,a^Yfj(a,a^) = fj(a,a j; )fj(a,a^ = 1. We note that 
{To is not defined on the whole space, which changes a domain and a range of U. 
Here we restrict each fj(a,a*) to one satisfying a relation 



f j (a,a^Nf j (a,J)=g j (N) 



for some function gj. For example, fj(a,a> 



y/(N+j-l)(N+j-2)-(N+l) 



i? . In this 



case 



D^D = U%DU 



because fj(a,a J< ) and the number operator N don't commute in general. 
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Let us show this with an example. For Bi in the two atoms case we consider a very simple 



case 



then it is easy to see 



where 



and 



U 



U 



i 



yW+T 



1 



v^V+T 



7WTi a J 



B x = UDU ] 



( x - a a 



^2(2^+3) VW+1^2(2N+3) ^2(2^+3) 



1 
V2 







1 

~V2 



2(2iv + r 



V y/2(2N-l) 

( 

D = 

V 

Compare this with (fHHJ) . Here we note that 



J_ a t _ vW^T Q t . 1 fl t 



Ny/2(2N-1) y / 2(2N-l) j 

\ 





2(2Ar + i) y 



In the two expressions Bi = UDU^ = UDW, each domain and range of U and U is different. 

The diagonal part D of B changes according to Uo, which is unavoidable due to the non- 
commutativity of operators in quantum physics. We call this phenomenon a U(l) ambiguity. 



5 Discussion 

We introduced the quantum diagonalization method and applied it to the (quantum) matrix B 
in (fT^j) (or A in Q) emerging from the Tavis-Cummings model and (re)obtained the explicit 
form of evolution operator for the one, two and three atoms case. To get the general case is 
not easy because of some technical reasons (numerical techniques are of course applicable). 
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Therefore, there are many applications to quantum optics or mathematical physics, see 
for example [4 . We can also apply the result to a quantum computation based on atoms of 
laser-cooled and trapped linearly in a cavity, see jlflj. 

We also make a comment on an application to a noncommutative differential geometry. 
From (}2"%)) we have a (quantum) unitary matrix U which gives the Maurer-Cartan forms 

L v = U^dU, R v = dUU-\ 

where d is some differential with respect to a and a* . These are fundamental objects in non- 
commutative chiral models. For the case of one, two and three atoms we can calculate the 
Maurer-Cartan forms exactly. Such a study is however beyond our scope of this paper. We 
expect that some researchers will develop the subject. 

We conclude this paper by making a comment. The Tavis-Cummings model is based on 
(only) two energy levels of atoms. However, an atom has in general infinitely many energy 
levels, so it is natural to use this possibility. We are also studying a quantum computation 
based on multi-level systems of atoms (a qudit theory) Therefore we would like to extend 
the Tavis-Cummings model based on two-levels to a model based on multi-levels. This is a 
very challenging task. 

Acknowledgment. We wish to thank Shin'ichi Nojiri for his helpful comments and suggestions. 



Appendix 

Proof of (gU) 
Here we show flUJ). The (l,l)-entry of Ue~ itgD U ] is 

+ e-^Va + e^"/ 4 u? 3 + ^'\ 2 U = 2u 2 n cos (tg^j + 2u? 2 cos (tg^j . 

because uf 3 = u\ 2 and u\ A = u\ x . Let us calculate both wf 1 and u\ 2 - 

2 {v^l-Z^ + v^ + T)} 2 
Un ~ 8(l + /5 2 )(l+ 7 2 ) 
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-^2 + 4 



1 f n , A V3{c(x - y) + b(x + y)} + ^ - ^ + 2bc 



8 1"'" /i 2 -z/ 2 



1 J -8N - 28 + 2^16(JV + 2) 2 + 9 
16(iV + 2) 2 + 9 



2iV - 7 + ^/16(iV + 2) 2 + 9 w+ (iV + 2) 



2^/16(A^ + 2)2 + 9 2^/d(iV + 2) 
where we have used the relations 

1 ji — u 1 x — y 1 fi + u 1 x + y 

P + ^ = —[— > 75~^ = ^ — ' 7 + - = , -- 7 = , 

pop o 707 c 

and (f4T)j) and (|41jl . Similarly, we have 



2 _ 2^ + 7+^16(^ + 2)2 + 9 _ v_{N + 2) 
-"11 _ 



2^/16(A^ + 2)2 + 9 2^/d(N + 2 

so that 



2u? 1 cos (^|) + 2w^ 2 cos f^-) 



U+ (iV + 2) _ / r — ; ~\ u_(iV + 2) 



cos tgJX+(N + 2) 7 cos WA_(iV + 2) = / 2 (iV + 2). 



2y/d(N + 2) V v / 2y/d(N + 2) 

The remaining 9 entries become more complicated because they contain a and (V . See [K] in 
detail. 
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